Stability Analysis of The Twisted Superconducting Semilocal Strings 

Julien Garaud and Mikhail S. Volkov 

Laboratoire de Mathematiques et Physique Theorique CNRS-UMR 6083, 
Universite de Tours, Pare de Grandmont, 37200 Tours, FRANCE 

We study the stability properties of the twisted vortex solutions in the semilocal Abelian 
Higgs model with a global SU(2) invariance. This model can be viewed as the Weinberg- 
Salam theory in the limit where the non-Abelian gauge field decouples, or as a two compo- 
nent Ginzburg-Landau theory. The twisted vortices are characterized by a constant global 
current Z, and for X — > they reduce to the semilocal strings, that is to the Abrikosov- 
Nielsen-Olesen vortices embedded into the semilocal model. Solutions withX / are more 
complex and, in particular, they are less energetic than the semilocal strings, which makes 
one hope that they could have better stability properties. We consider the generic field fluc- 
tuations around the twisted vortex within the linear perturbation theory and apply the Jacobi 
criterion to test the existence of the negative modes in the spectrum of the fluctuation op- 
erator. We find that twisted vortices do not have the homogeneous instability known for 
the semilocal strings, neither do they have inhomogeneous instabilities whose wavelength is 
less than a certain critical value. This implies that short enough vortex pieces are perturba- 
tively stable and suggests that small vortex loops could perhaps be stable as well. For longer 
wavelength perturbations there is exactly one negative mode in the spectrum whose growth 
entails a segmentation of the uniform vortex into a non-uniform, 'sausage like' structure. 
This instability is qualitatively similar to the hydrodynamical Plateau-Rayleigh instability 
of a water jet or to the Gregory-Laflamme instability of black strings in the theory of gravity 
in higher dimensions. 
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I. INTRODUCTION 

Vortices have been the subject of intense studies ever since their discovery by Abrikosov in 
the context of the Ginzburg-Landau theory of superconductivity and by Nielsen and Olesen in 
relativistic field theory QJJJ - They find numerous applications in many domains of physics ranging 
from high energy physics and cosmology Q2Q to various branches of condensed matter physics, 
such as superconductivity [3] and superfluidity |4|] models. In these applications the vortices are 
most often considered within the original model of Abrikosov, Nielsen and Olesen (ANO) yj], 
which means as solutions of the Abelian Higgs model containing an Abelian vector coupled to 
a complex scalar $. However, vortices can also exist in others, more general field theory models. 

A particular field theory containing an Abelian vector and two complex scalars $1 and $ 2 
with a global SU(2) invariance, dubbed semilocal model, has been extensively discussed recently 
I5Q- This model can be considered as a 'minimal' generalization of the original ANO theory with 
the following interesting properties. First, it can be viewed as the special limit of the electroweak 
theory of Weinberg and Salam in which the weak mixing angle is n/2 and the non- Abelian gauge 
field decouples (a). Secondly, it describes states with a multicomponent order parameter in con- 
densed matter physics, as for example it the two band superconductivity models or in superfluid 
systems [6]. Finally, it can be regarded as sector of supersymmetric field theories [71] or as part 
of one of the Grand Unification Theories that could presumably describe cosmic strings in the 
early Universe [8D. All this implies that the semilocal model is likely to have interesting physical 
application, which has inspired considerable interest towards this theory and its solutions. 

The original ANO vortex can be embedded into the semilocal theory, since identifying its field 
A^ with that of the semilocal model and setting also $i = $ and $2 = solves the field equations 
of the theory [9]. Such an embedded vortex is called semilocal string and it behaves identically to 
the original vortex as long as <3> 2 = 0, but its properties can become quite different as soon as the 
scalar $ 2 is excited [5]. For example, although the original ANO vortex is stable, fluctuations of 
$ 2 render its embedded version unstable in the parameter region where the Higgs boson mass is 

o 

larger than the vector boson mass [10]. 

Not all vortices in the semilocal theory should necessarily be of the embedded ANO type. It has 
been known for some time that more general vortex solutions (sometimes called 'skyrmions') exist 
in the theory in the special limit where the Higgs boson mass is equal to the vector boson mass 
Jlp| . [U]. For these solutions $i and $ 2 are both non-vanishing. Quite recently vortex solutions 
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with a similar property have been constructed within the whole parameter region of the theory in 
which the Higgs boson mass is larger than the vector boson mass II 1211 . They are characterized by 
the twist: a ^-dependent relative phase exp(iKz) between $1 and $ 2 . The twist gives rise to a 
constant global current along the vortex, whose value, X, is a parameter of the solutions. For this 
reason the solutions are called twisted superconducting semilocal strings, or twisted vortices for 
short. For X — > they reduce to the semilocal strings but for X ^ they are physically different. 
In particular, it turns out that solutions with X ^ are less energetic as compared to the semilocal 
strings. This suggests that, since they are energetically favoured, it is predominantly the twisted 
vortices and not the semilocal strings which should be created in physical processes leading to 
vortex formation. Therefore, as far as physical applications is concerned, twisted vortices could 
be more important than the semilocal strings. However, since they exist precisely in the same 
parameter region where the latter are unstable, the twisted vortices may well be also unstable, 
which would somewhat delimit their importance. On the other hand, one may argue that they 
should have better stability properties than the semilocal strings, since they are less energetic. 

Motivated by this, we study in this paper the stability of the twisted vortices within the linear 
perturbation theory by analysing the linearized equations for fluctuations around the twisted vor- 
tex background. Decomposing the perturbations into a sum over Fourier modes proportional to 
exp{i(ujt + Kz + rmp)}, the variables in the fluctuation equations separate and the equations reduce 
to a Schrodinger type spectral problem with the eigenvalue cu 2 . If this problem admits eigenstates 
with cu 2 < then the perturbations will be growing in time and so the background will be unstable. 
In order to find out whether such negative modes exist, we apply the Jacobi criterion [[130, which 
only requires the knowledge of the u = solutions of the fluctuation equations. 

We arrive at the following conclusions. For 1^0 the perturbation equations do not admit 
negative mode solutions independent of z. This means that the twisted vortices do not have the 
homogeneous instability known for the semilocal strings IllOn . Moreover, considering inhomoge- 
neous 2-dependent perturbations of the fundamental twisted vortex, we find that the corresponding 
eigenvalues are always positive, apart from those corresponding to the modes with m = and with 
the wavelength 

7T 

A > A m j n = n^p (i) 

where K is the twist parameter of the vortex. It follows then that short vortex pieces obtained 
by imposing the periodic boundary conditions on the infinite vortex will be perturbatively stable 
if their length, L, is less then A min . If L > A min then the vortex will have enough room to ac- 
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commodate inhomogeneous instability modes whose growth will lead to its fragmentation into a 
non-uniform, 'sausage like' structure characterized by zones of charge accumulation and by an 
inhomogeneous current density. One needs to go beyond the linearized approximation in order to 
fully understand the development of this instability. For higher winding number vortices we find 
additional instabilities corresponding to their splitting into vortices of lower winding number. 

Our general conclusion is that the twisted vortices indeed have better stability properties than 
the semilocal strings since, unlike the latter, they can be stabilized by imposing periodic boundary 
conditions. This suggests that small stationary vortex loops, if exist, might be stable, which opens 
intriguing research perspectives that may lead to interesting applications in various domains of 
physics ranging form condensed matter physics to cosmology. 

The rest of the paper is organized as follows. In Sec.II we summarize the essential properties 
of the twisted vortices. Sec. Ill contains the analysis of generic perturbations around the twisted 
vortex background in the linearized approximation: the mode decomposition, separation of vari- 
ables in the perturbation equations, gauge fixing, and the reduction to a Schrodinger type spectral 
problem. The existence of the negative modes in the spectrum is demonstrated with the use of the 
Jacobi criterion in Sec.IV, where these modes are also explicitly constructed and their dispersion 
relation is determined. The physical manifestations of the instability are discussed in Sec.V, where 
analogies from other branches of physics are also discussed. We summarize our results in Sec. VI, 
while the complete system of perturbation equations is presented in the Appendix. 

The issue of vortex perturbations is not so often considered in the literature. Although the 
basic stability pattern for the ANO vortices was understood long ago [14], the systematic stability 
analysis in this case was carried out only relatively recently jl5\\. In our analysis we essentially 
follow the approach of IU5I1 . although working in a different gauge. 



II. THE TWISTED SEMILOCAL VORTICES 



A. The semilocal model 



The model considered here is obtained by replacing the complex scalar of the usual Abelian- 
Higgs model by a doublet of complex scalars, $ tr = ($1, $2)- This theory, the semilocal model, 
is by construction invariant under the internal symmetry SU(2) g i ob ai x U(l)i oca i- It is worth noting 
that this semilocal model describes the bosonic sector of the electroweak theory in the limit where 
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the weak mixing angles is n/2 and the non-Abelian gauge field decouples 

The Lagrangian density in terms of suitably rescaled dimensionless coordinates and fields reads 

C = -\f, u F^ + (D^) ] - | - l) 2 , (2) 

where F^ u = d^A v — d u A^ and D^<& = d^<& — iv4 M 3>, the spacetime metric signature being 
(+, — , — , — ). The spectrum of this theory consists of a massive vector boson whose mass in the 
dimensionless units chosen is m v = a/2, of two massless Goldstone scalars, and of a Higgs boson 
with the mass m H = \fj3 m v . The fields transform under SU(2) g i ob ai x U(l)i oca i as 

A^A^ + d,A(x), $^U$, (3) 

with U = e lA ( x )+^aT a w h ere Ta are the Pauli matrices and 9 a are constant parameters. The Euler- 
Lagrange equations of motion obtained by varying the Lagrangian ((2|) read 

d^ = i{(A,$) f 

= -15 ($ f $ - l) $, (4) 

where the dagger stands for the hermitian conjugation. 

The internal symmetries of the theory give rise to several independently conserved Noether 
currents, one of which will be important in what follows, 

J M = U(m* 2 D^ 2 ). (5) 



B. The twisted superconducting vortices 

Eqs.© admit stationary, cylindrically symmetric solutions of vortex type II 1211 . For these solu- 
tions the fields are parametrized in cylindrical coordinates as 

I fi(py Ntp \ 

A.dx^ = a 2 (p)(Qdt + Kdz) + NaAp^dkp, , (6) 

I fJp) e i{Vt+Mip+Kz) J 

where all functions of p are real and N, M are two integer winding numbers. The two real pa- 
rameters and K can be seen, respectively, as the relative rotation and twist between the two 
components of the scalar doublet. With this parametrization Eqs.© reduce to a set of four non- 
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linear ordinary differential equations, 
1 



P 



(pa' 2 )' = 2a 2 f 2 + 2(a 2 -l)fl 



~ N ) 



P (si\ = 2/l 2 ( ai - 1) + 2f 2 L 

- p (pf[Y = h { ^ 2{ai ~ + (K 2 - tf)a\ + p (A 2 + f! - l)) , 

- p M)' = h ( (M ~f Ql)2 + (K 2 - fi 2 )(a 2 - I) 2 + /? (fl + /f - l)) , (7) 

with ' = ■£-. One can consistently set in these equations / 2 = a 2 = 0. The remaining two 
non-trivial equations for / l5 a x reduce then to the ANO system [1] and so the solutions will be 
the ANO vortices embedded into the semilocal theory. Such embedded solutions are sometimes 
called semilocal strings |5|]. 

In addition, for @ > 1, Eqs.fT]) possess also more general solutions, called twisted vortices, 
which have f 2 ^ and a 2 ^ ftl2Q. Their numerical profiles (see FigHJ) show that / a , ai behave 
qualitatively in the same way as for the ANO vortex, while f 2 , a 2 develop non-zero condensate 
values in the vortex core and tend to zero for p — > oo. For a given (3 > 1 these solutions comprise 
a three parameter family labeled by N = l,2,...,byM = 0,l,...iV — 1, and by a real parameter 
q = --^j/r | M (0): the value of the M-th derivative of f 2 at the vortex center. In what follows we 
shall call q condensate parameter and shall be mainly considering solutions with M = 0, so that 
Q = .MO). These parameters determine the value of the combination K 2 — f2 2 which turns out 
to be positive for all twisted solutions. It is worth noting that the ansatz (J6]) preserves its structure 
under Lorentz boosts along the vortex axis - if we assume that (O, K) transform as components 
of a spacetime vector. Since the Lorentz invariant norm of this vector, K 2 — Vt 2 , is positive, the 
vector is spacelike, and so one can boost away its temporal component by passing to the restframe 
where Vt = 0. On the other hand, the twist K is an essential parameter that cannot be removed, 
which is why the vortices are called twisted. Since the equations are invariant under K — > —K, in 
what follows we can assume without loss of generality that K > 0. 

Associated to the twist there is a physical parameter: the total current © through the vortex 
cross section, 

1= [ J z d 2 x= [ d 2 x(A z - K)^ 2 . (8) 



This is zero for the embedded ANO solutions and non-zero for the twisted vortices. The current 
is a function of the condensate parameter q (also of N, M). In the limit q — >• one has f 2 — > 0, 



Figure 1: Left: profile functions for the twisted vortex solution with (3 = 2, N = 1 and q = 0.5. Right: 
the restframe (Q = 0) vortex energy E, current T, and twist K against the condensate parameter q for the 
N = 2, M = twisted solutions with (3 = 2. 



a 2 — > 0, the current vanishes and the twisted vortices reduce to the semilocal strings. For q ^ the 
current is non-zero and it becomes arbitrarily large when q — > 1 (although this is not so visible 
in FigCQ). Interestingly, both the twist K and the vortex energy per unit length, E = J T®d 2 x, 
decrease when the current increases (see Fig{l]), for large currents the energy approaching the 
lower bound 27r|iV| lll2n . The twisted current carrying vortices are thus energetically favored 
when compared to the currentless semilocal strings, and so that one may expect that they could 
have better stability properties than the latter. 



IH. PERTURBATIONS OF VORTICES 



In order to investigate the stability of twisted vortices, one has to examine the dynamics of 
their perturbations. Let \I/(f) collectively denote the fields of a background vortex solution. Let us 
consider small perturbations around this background, 

^(r) -> ^(r) + 5^(r)e iujt . 

Inserting this into the field equations © and linearizing with respect to the perturbations, one can 
put the linearized equations into a Schrodinger type eigenvalue problem form, 

(-A + U)6V = uj 2 5V, (9) 
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where the potential U is determined by the background fields ^ . If the spectrum of the differential 
operator in the left hand side of this equation is positive, then all the eigenfrequencies u are real 
and so that the background ^ is linearly stable. If on the other hand there are bound states with 
to 2 < 0, the frequency uj will be imaginary and the corresponding perturbation modes will grow in 
time as e^' thus indicating the instability of the background. 



A. Special fluctuations around the ANO vortex 

Let us first briefly consider the q — > limit when the background corresponds to the embedded 
ANO vortex. It is easy to analyze a particular type of perturbations around this solution by simply 
choosing the amplitudes f 2 and a 2 in Eq.© to be small, such that the fields read 

(A^ + 5A^)dx^ = 5a 2 (p)(cudt + kdz) + Na x (p)dLp, 
( 



$ + 5$ 



(10) 



Here we have replaced fi, K, M, respectively, par cu, k, p in order to emphasize that these param- 
eters relate to the perturbation and not to the background. Of course, this gives only a special type 
of perturbations and not the most general one. The perturbation equations are obtained by simply 
linearizing Eqs.CZ]) with respect to f 2 = 5f 2 and a 2 = 5a 2 . The resulting linear equation for 5a 2 
decouples and one can show that its only bounded solution is 5a 2 = 0. The equation for 5f 2 reads 

--Wm' + [ {Nai ~ V? + P (fl - 1)) 5f 2 = e5f 2 , (11) 

where fi and a\ refer to the iV-th background ANO solution, 77 = 0, 1, . . . iV — 1, and 

£ = u 2 -k 2 . (12) 

The numerical analysis reveals flp| that if /3 > 1 and for = 1, 77 = there is a bound state 
solution of Eq. (fTTj) . 

6f 2 (p) = Mp), e = e (P)<0, (13) 

where ^(0) 7^ 0, ^0(00) = and J °° \ ^f \ 2 pdp = 1. The eigenvalue e is related to the restframe 
value of the background twist K in the q — ► limit (see FigQ]) as lll2ll 

e = -K 2 . (14) 
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The frequency therefore satisfies the dispersion relation 

u 2 = k 2 - K 2 (15) 



and so the perturbations 



will be growing in time if 



8<& 2 = e iwt+ikz ^ (p) (16) 



\k\ < \K\. (17) 



The ANO vortices, when embedded into the semilocal theory, are therefore dynamically unstable 
I loll . The modes (fT6l) depend on z and so they describe inhomogeneous instabilities, apart from 
the k = mode corresponding to the homogeneous instability of the embedded vortex II 1011 . For 



1*1 > \K\ (18) 

one has to 2 > and Eq. (fT6l) gives stationary deformations of the ANO vortex. Eq.dlOl) (with 5a 2 = 
0) then describes the ANO background plus the first order correction due to the twist/current. For 
k = ±K one has uj = and so Eq.(fT6l) gives zero modes corresponding to static restframe 
deformations of the ANO vortex by the current. 

B. Generic perturbations of twisted vortices 

Since the twisted vortices are unstable in the q — > limit, when their current vanishes, one 
can expect that they will presumably be unstable also for q ^ 0, at least for q <^ 1 when the 
current is small. However, it is not logically excluded that they may become stable for larger 
values of the current. Let us therefore consider small fluctuations around the generic twisted 
vortex configuration, 

$ = $M + s^°\ = Af + 8Af, (19) 

where <fr^, A^ are given by Eq.©. In order to analyze the dynamics of perturbations it is con- 
venient to make use of the fact that the background configurations are stationary and cylindrically 



symmetric 1121]. Since the corresponding symmetry generators d/dt, d/dz, d/dcp commute be- 
tween themselves, there exists a gauge where the background fields do not depend on t, z, (p. 
However, one cannot pass to this gauge using only the local U(l) gauge symmetry that we have 
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in our disposal, since adjusting one gauge function A(x) in Eq.® does not allow to eliminate the 
phases of the two Higgs field components at the same time. However, the following technical trick 
can be employed. 

Let us introduce an auxiliary SU(2) gauge field = r a W^ which is pure gauge. Replacing 
then in Eqs.©,© the covariant derivative as 

D, - V, = D, - iW«r\ (20) 

the SU(2) g i oba i x U(l)iocai symmetry © can be promoted to the SU(2)i oca i x U(l)i oca i gauge 
transformations: 

$^U$, A^ + W^ViA^ + W^ + id^V- 1 , (21) 
where U = e lA ( x )+^a(x)r a _ § mce [ s p Ure gauge, one can gauge it away and then we return 
to the previous formulation of the theory. However, allowing for non-zero values of gives 
us an additional local SU(2) gauge freedom, which can be useful. Although within the original 
semilocal model the field is merely a technical tool, and so we shall have to remove it at the 
end of calculations, it can be viewed as a physical field if the semilocal model is considered as the 
limit of the Weinberg-Salam theory. 

Let us now apply to Eq.(|T9T) the gauge transformation (I2TT) generated by 

-iNip o 

U = . (22) 

q e -i{Q.t+Kz+M<p) 

This transforms Af + 5Af -> A^ ] + SA 1 ? and $M + $W + as well as = 

-> where the new background fields 

*W = , W^da? = —(Qdt + Kdz + (M — N)d<p), 

\fMJ 2 

AWdbf = (a 2 (p) - l -)(Sldt + Kdz) + (N( ai (p) - I) - *i)dtp, (23) 
while the old and new perturbations are related as 

Mjl = 5Af, 5$S 1] = e- iN ^ ] , = e -Km+Kz+M^[o] (24) 

As a result, the background fields depend now only on p. The price we pay for this is a new field 
Wji that appears in the new gauge. Moreover, we observe that the azimuthal components of the 
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vector fields, and W}p , do not vanish at p — 0, which means that the fields are not defined at 
the symmetry axis. Nevertheless, we can work in this gauge and calculate 8Ap, 5$$, provided 
that when transformed back to the old gauge where the background fields are globally regular, the 
perturbations 5A$ and obtained via Eq. (l24j) will also be regular. 

Inserting $ = 3>M + <5<frM an d A^ = A^ + 5A$ to Eqs.©, replacing the covariant derivatives 
by Vp, linearizing with respect to the perturbations and omitting the superscript gives 

d^8A u - d u d,5A^ = 

= i [(D„*)t(?$ + (V v 6&)*& - - 6&V U &] - 2\<$>\ 2 5A l/ , 

V^WSQ - 2i5A^<f> - i^d^A^ = -f3 (2|$| 2 - l) 5$ - /M$ f $ 2 , (25) 

where P M = — i(A^ + W^) and the background fields <fr, A^, are given by (|23T) . 
These equations are invariant under the U(l) gauge transformations 

SAp^SAp + dnSkix), 5® ^ 5$ + i$5A(x) } (26) 

which is the infinitesimal version of the local U(l) gauge symmetry contained in (12TT) . 

C. Separation of variables 

Since the coefficients in Eqs. (|25l) depend only on p, it is straightforward to separate the variables 
in these equations by making the Fourier type mode decompositions, 

5$ a = cos M + mi f + Kz ) (C A,m (p) + iC' K ' m (p)) 

+ sm(ut + mp + kz) K' K ' m (p) + i^' K ' m (p)) , 
SAfj,, = &' K ' m (p) cos M + rrup + kz) + X ^' m (p) sin(u;t + rmp + kz), (27) 

Lj,K,m 

where a = 1, 2 and p = 1, 2, 3, 4 with x^ = (t, p, z, cp), respectively. One can similarly decompose 
the gauge function SA in (|26l ), 

5A= (cos(ut + rmp + Kz)a"' K ' m {p) + sin(ut + rmp + Kz)^' K ' m (p)), (28) 
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and so that the gauge transformations (|26|) assume the form 



u),K,m 
a 





! 


w,K,m 
"a 


' J a i i 


s-uj,K,m, 
SI 




X2 




£LO,K,m 

S3 




S4 





V 




ra 


Ya 1 J a" 


Xi - 




S2 ~~ 


+ £2 + (" ) 




u>,K,m u},K,m 
A3 ~~ Ka 


uj,K,m 

Xi - 


-> X4 A,m - mc^' K ' m 



SI ' 

X 2 

j-w,K,,m 
S3 

g^m ^ ^,„, m + ^y^m ^ _ maWl , h m | (29) 

where ' = 4-. 

dp 

Inserting decompositions (|27l) to the equations (1251) . the t, z, cp variables separate and for fixed 
values of u>, k, m one obtains a system of 16 ODE's for the 16 real functions 0a ,ft,m (p), vr^ ' K ' m (p), 
u a' K ' m (p)> ' l l ) a' K '' m {p)'> Xp' K ' m (p)i 4>' K ' m (p)- A further inspection reveals that these 16 equations 
actually split into two independent subsystems of the same size. Specifically, the 8 amplitudes in 
the left column in (l29~l) . whose transformations involve only the gauge functions 7 w ' K > m , satisfy a 
closed system of 8 equations. Similarly, the amplitudes in the right column in (l29l) satisfy a closed 
system of 8 equations. These two groups of equations become identical upon the replacement 

n a ~^ 'Pa i Va ~ > V a 1 S2 ~^ ~ A2 ) 

x;' K ' m ^C, K ' m , (P = 1,3,4). (30) 

As a result, without any loss of generality we can restrict our consideration to the 8 equations 
of the first group. They are explicitly listed in the Appendix, Eqs. dA. II )— (1A.8I) . where we have 
omitted the superscripts to, k, m to simplify the notation. As we shall see later, the two groups of 
equations describe simply the real and imaginary parts of the perturbations. 



D. Gauge fixing and reduction to a Schrodinger form 

One can check that Eqs. (|A. II )— (IA.8I) are invariant under the gauge transformations expressed by 
(l29~l) . In addition, they satisfy the differential identity (|A.9I) expressing the fact that the divergence 
of the right hand side of the SA^ equation in Eqs.(l25l) should vanish, since the divergence of its 
left hand side vanishes identically. The identity (|A.9I) expresses one of the 8 equations in terms of 
the remaining ones, and so only 7 equations are actually independent. As a result, one can exclude 
one of the equations from consideration. 
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Let us now specialize to the case of purely magnetic backgrounds by setting 

Q = 0. (31) 

There is no loss of generality whatsoever in imposing this condition. Indeed, since the twisted 
vortices with £7 ^ can be obtained from those with = by a Lorentz boost, the same should 
apply to their perturbations. It is therefore sufficient to consider perturbations in the f2 = case, 
since Lorentz boosting them will give perturbations for the ^ backgrounds. 

Next, we fix the gauge freedom (I26l),(l2~9l) by imposing the temporal gauge condition 

<L4 Q = ^ = ^' K ' m = 0. (32) 

This fixes the gauge completely if uj ^ 0, although leaving a residual gauge freedom in the uj = 
sector expressed by (T26l) with 5 A independent on time. Now, equation (IA.7I) for £x is the Gauss 
constraint, and when ^ = it assumes the form 

(I 77? \ 

X2 + -X2-^U- <3 - 2f x u x - 2f 2 u 2 j = 0. (33) 

We use this equation to express £3 in terms of the other variables, 

6 = - U + ^ - ^ - 2(M + m) , 04) 
« V P P J 

which is possible as long as k 7^ 0, and so that from now on we can exclude £ 3 from the consid- 
eration. Since one of the 8 equations is redundant, we can also exclude from the consideration 
the equation (IA.5I) for £ 3 . As a result, we now have only 6 equations for the 6 independent field 
amplitudes <p a , v a , % 2 , ^4 (see the Appendix for more details). 

Remarkably, the elimination of £3 also truncates the residual gauge freedom of time indepen- 
dent gauge transformations. Specifically, the gauge transformation (|29l) with u = leave all the 
equations invariant, in particular the Gauss constraint (|33l) . However, the gauge invariance of the 
latter is only insured by the overall factor of to, expressing the fact that Eq. (|33l) is a total time 
derivative. When uj = 0, the whole expression in (1331) is zero and so that it is gauge invariant. 



However, the expression between the parenthesis in (1331) . sometimes called 'strong Gauss con- 
straint' lllql . is not gauge invariant under (|29| ) for an arbitrary gauge function 7 = 7°' K > m but only 



if 7 satisfies the condition 



7 ,/ + -7 / =^ + 2(A 2 + /|) + k 2 )7. (35) 



P 
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Since we have used the strong Gauss constraint to express £ 3 , we have actually reduced the residual 
gauge freedom to the two parameter family of solutions of this differential equation. It is not then 
difficult the see that if 7 does not vanish identically, then it must be unbounded for small or for 
large p (or in both limits). The behavior of the pure gauge modes produced by this 7 is thus 
incompatible with the boundary conditions at the origin or at infinity (see Eq.(f38b.(l39l) below), 
unless 7 = 0. The only admissible solution of Eq.d35l) is therefore 7 = 0, which fixes the residual 
gauge freedom completely. 

We finally arrive at a system of 6 independent second order equations. They are listed in the 
Appendix, Eqs. dA. 141) — (IA. 191) . and they can be rewritten in the form 

- + = u 2 ^! . (36) 

Here \& is a 6-component vector whose components are expressed in terms of a , u a , % 2 , £ 4 and 
their first derivatives, and U is a potential energy matrix determined by the background solution 
and depending also on k, m. This determines a linear eigenvalue problem on a half-line, p E 
[0, 00), the points p = 0, 00 being singular points of the differential equations. 



E. Boundary conditions for perturbations 

Let us study the local behavior of solutions of Eqs. ((36l) for small and large p. It is convenient 
to introduce the linear combinations 



"a ! 



<p a = K + h a , v a = h+ - h a 

U = p(g + -g~), X2 = g + + g~- (37) 



Using the asymptotic expansions of the background solutions at small p [|12|1 to determine the 
asymptotic behavior of U one can construct series solutions of Eqs. (|36l) in the vicinity of p = 0. 
One finds that bounded at p = solutions behave as 

hf = Af p lNTml + . . . , hf = Afp\ M ^ + . . . , g ± = B ± p\ m ^ 1 \ + . . . , (38) 

where A^, B ± are integration constants and the dots denote the subleading terms. 

Let us now consider the asymptotic region, p — > oo. Setting the background field amplitudes 
to their vacuum values, f 1 = a± = 1, f 2 = a 2 = 0, so that the background fields (1231) become 
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pure gauge, Eqs. dA. 14b — (IA. 191) decouple from each other and read (with T> + being defined in the 
Appendix) 

(tti 2 \ 
-T>+ + + K 2 - uj 2 + 2(3 \ 0! = 0, 

( 777^ \ 
-V + + — + k 2 -uo 2 + 2)^ = 0, 

(-A + {N - M / m) \ ( K ±KT-^hi = o. 

The first equation here represents a massive Higgs boson with m H = \/2/?, the next three corre- 
spond to one longitudinal and two transverse degrees of polarization of a massive vector boson 
with m v = a/2, while the last two equations for describe a pair of Goldstone particles. The 
fluctuation equations therefore reproduce correctly the physical spectrum of the semilocal model, 
and the asymptotic behaviour of the solutions for p — ► oo is 

0i = exp{-/i H p} + . . . , !/! = _! exp{-/i„p} + . . . , 

9 ± = -p: exp{-p v p} + • • • , = — exp{-p ± p} + . . . , (39) 

Vp Vp 

where ai,a 2 ,a±, b± are integration constants and 

pl = K 2 -uj 2 + 2(3, P 2 v = k 2 -uj 2 + 2, p 2 ± = (k ± K) 2 - uj 2 . (40) 



IV. STABILITY TEST 

Summarizing the above analysis, we have arrived at the eigenvalue problem (1361) . and now we 
wish to find out whether it admits bound state solutions with uj 2 < 0. If exist, such solutions would 
correspond to unstable modes of the background vortex configuration. 

Bound state solutions of Eqs.(l36l) should be everywhere regular and so they satisfy the boundary 
conditions (|38l) ,d39l) at the origin and at infinity. One possibility to proceed would then be to 
directly integrate Eqs.d36l) looking for solutions with uj 2 < 0, which would require solving the 
boundary value problem for the 6 coupled second order differential equations. Fortunately, there 
is a faster way, since to know whether the negative modes exist or not it is not actually necessary to 
construct them explicitly. A simple method that we can apply to reveal their existence is to use the 
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Jacobi criterion [13] (see II 1711 for applications of this method for the monopole stability problem). 
This method essentially relates to the well known fact that the ground state wave function does 
not oscillate, the first excited state has one node and so on. It follows then that if the zero energy 
solution of the Schrodinger equation oscillates, then the ground state energy eigenvalue is negative. 
When applied to multichannel systems as in our case, the Jacobi method gives the following recipe. 

Let \I/( s )(p) where s = 1, ... 6 be 6 linearly independent and regular at the origin solutions of 
Eqs. (|36l) with cu 2 = 0, each of them being a 6-component vector ifflJp). These solutions can be 
chosen to satisfy, for example, the conditions ty 1 ,^ (p ) = $1, where p is a point close to the origin. 



Equally, each of them can be obtained by taking the boundary conditions (1381) . setting to zero 5 out 
of the 6 integration constants in (|38l) , and integrating numerically Eqs.(|36l) with uo 2 = towards 



large values of p. Now, if the determinant 



A(p) 



*0)(P) ■ ■ *(6)(P) 



(6)' 



(6)V 



(41) 



vanishes somewhere, then there exists a negative part of the spectrum. Thus we can decide whether 
the backgrounds is stable or not by simply calculating the determinant (|4TT) - a much easier task 
than solving the boundary value problem (|36l) to directly find the eigenvalues. According to [18], 
the number of instabilities is equal to the number of nodes of A(p). 



A. Jacobi test in the ANO limit 



Before studying the general case, let us consider the limit where q = f 2 = a 2 = and the 
twisted vortices reduce to the embedded ANO solutions. We already know that in this case there 
is a particular perturbation (fTOl ) leading to the eigenvalue problem (fTTT) admitting the negative 
mode (fT6l) . Now we shall be able to take into account all the remaining perturbation modes. The 
perturbation equations in the ANO limit split into the system of four coupled equations (IA.21I) 
plus two decoupled equations (|A.22I) . 

The four equations (IA.21I ) do not contain perturbations of the second component of the Higgs 
field and so they actually describe the dynamics of the perturbed ANO vortex within the original 
one component ANO model. We therefore expect that for N = 1 these equations do not admit 
bound state solutions with lo 2 < 0, since the ANO vortex is stable within the original ANO theory. 
However, since for (3 > 1 the multivortex configurations of higher winding numbers are unstable 
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Figure 2: Jacobi determinant Ajv ;m for solutions of Eqs. (IA.21l) with lo = n = and h\ given by (IA.22I ) 
with lo = m = 0, k = K in the case of the = 2 ANO background with N = 1. Since Ai m do not 
vanish, the iV = 1 ANO vortex is stable within the original ANO theory. The vanishing of A2,2 indicates 
the splitting instability of the N = 2 vortex. The vanishing of h\ shows the instability of the N = 1 vortex 
embedded into the semilocal theory. 



with respect to splitting into vortices of lower winding numbers Ill4ll . lll5h . Eqs. (IA.21l) should admit 
bound state solutions with uj 2 < for j3 > 1 and for N > 2. We can now test these assertions using 
the Jacobi criterion. Integrating Eqs. (IA.21l) with uj 2 = with the regular at the origin boundary 
conditions (|38l) to construct their four linearly independent solutions, we calculate the determinant 
A(p) as in Eq.(l41~l). It turns out then that A(p) is indeed always positive for N — 1, while for 
N = 2 and j3 > 1 it changes sign for quadrupole perturbation modes with m = 2 (see Fig©. All 
this agrees with the expected results. 

Let us now consider the two decoupled equations (|A.22I) . They have exactly the same structure 
as Eq.(fTTT). up to the replacement 



Sf2 



7] — >• M ± m, k — > K ± k 



(42) 



We therefore immediately know that for iV = 1, in which case M = 0, choosing m — 0, 
these equations admit bound state solutions. We can also directly apply the Jacobi criterion to 
Eqs. (IA.22l) . Setting uj = and also K + k = (or K — n = 0) and integrating gives a solution 
for h\ (or for h^) that passes through zero once (see Fig[2]) thus showing that there is exactly one 
bound state. 



18 



B. The N = 1 twisted vortices 

We now have all the necessary ingredients to check the stability of generic twisted vortices. To 
this end we integrate numerically the system of 6 coupled equations (1A. 141) — (IA. 191) with to 2 = 
with the regular at the origin boundary conditions d38l to construct the Jacobi determinant. We 
consider first of all the case of the fundamental N = 1 twisted vortex, and for the perturbation 
parameters we choose the same values which give rise to the instability in the ANO limit: m = 
and k = K. The resulting Jacobi determinant A(p) for several values of the background conden- 
sate parameter q is shown in FigfH We find that not only for q — > (the ANO limit) but also for 
all other values in the interval < q < 1 the determinant crosses zero. All these solutions are 
therefore unstable with respect to axially symmetric perturbations, and we have gone up to (5 = 10 
to check that the number of instabilities, given by number of zeroes of A(p), is exactly one. 




Figure 3: Right: Ai ( q for k = K for the twisted backgrounds with (3 = 2 and q G [0; 0.35]. The vanishing 
of Ai o indicates the existence of a negative mode. Left: perturbations of the N = 1 twisted vortex in the 
m > sector. No instabilities appear. Here (3 = 2, q = 0.5, k = K. 

We have also checked the sectors with m > but found no instabilities there. In Figj3] the 
Jacobi determinants for m = 1, 2 are shown, they are everywhere positive. This is not surprising, 
since increasing m increases the centrifugal energy thus rendering less probable the existence of 
bound states. Summarizing, all instabilities of the fundamental iV = 1 twisted vortex reside in the 
axially symmetric m = sector. 
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C. Solutions with N > 1 

Vortices with N > 1 have additional instabilities. First, they have the same instability in the 
m = sector as the N = 1 solutions. This can be checked by calculating the Jacobi determinant 
for the m = perturbations; see FiglH In addition, solutions with N > 1 can be unstable with 
respect to splitting to vortices of lower winding number. 

The latter instability is present already in the ANO limit (see Fig©. In the ANO case the 
existence of this splitting instability can also be inferred from the energy considerations: in the 
simplest case is suffices to compare the energy of the N = 2 vortex to the doubled energy of the 
N — 1 vortex. It turns out then that the former is higher than the latter for f3 > 1, which means 
that the vortex splitting is energetically favoured. 




Figure 4: Left: perturbations of the N = 2 vortices. Here (3 = 2, k = K and q G [0;0.35]. The same 
m = instability as for the N = 1 vortex is detected. N = 2 vortices are however stable with respect to 
perturbations with m > 0, m / 2. Here [3 = 2, k = K and q = 0.5. Right: they are unstable with respect 
to the m = 2 quadrupole deformations - splitting instability. Here j3 = 2, k = K/2 and q G [0; 0.35]. 

A similar energy argument can also be applied to the twisted vortices. However, the consider- 
ation in this case is complicated by the fact that, apart from the winding number N, the twisted 
vortices also carry additional parameters, in particular the current X. One can also compare the 
energy of a N = 2 twisted vortex, E(N = 2,1), with a sum of energies of two N = 1 vortices, 
E(N = l,2i) + E(N = 1,X 2 ), but then one should consider different possibilities for values of 
X\ and X 2 . Assuming that X\ + X 2 = X and choosing X\ = X 2 = X/2 one finds that the splitting 
is indeed energetically favoured if X is small, but not for large enough X II 1211 . However, one can 
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also consider the case where X\ is large and positive, while X 2 is large and negative, their sum 
Xi + X 2 = X being fixed. Since the energy of twisted vortices decreases with the current (see 
FiglH), it follows that whatever the value of X is, one can always adjust Xi and X 2 such that the 
splitting will be energetically favorable. 

We therefore expect all twisted vortices with N > 1 to be unstable with respect to the splitting. 
This expectation is confirmed by calculating the Jacobi determinant in the m = 2 sector, since for 
all N = 2 backgrounds we tested the determinant vanishes somewhere; see FigJH We have also 
checked that N = 3 vortices exhibit the splitting instability for m = 2,3. 

Summarizing what has been said, N > 1 current carrying vortices exhibit the same axially 
symmetric m = instability as the fundamental N — 1 vortex. In addition they are unstable with 
respect to splitting into vortices of lower winding number. We have not found negative modes in 
sectors with m = 1 or for m > N. 

D. Explicit computation of the eigenvalue 

Having revealed the existence of the negative modes in the spectrum of the fluctuation operator, 
we have also managed to construct them explicitly. Such a construction is considerably more 
involved than applying the Jacobi criterion, since it requires solving the boundary value problem 
for the 6 coupled equations (1 A . 1 41) — (I A .191) with the boundary conditions ((381) and d39l) . In addition, 
one should solve at the same time the 4 equations © to generate the background profiles. 

It turns out that for the N = 1 background Eqs. dA. 141) — (|A. 19b admit two different bound state 
solutions in the m = sector with eigenvalues u?^ and u%y where uj^Jk) < for k E {—2K, 0) 
and u>? 2 Jk) < for k E (0, 2K), also uZJk) = uj 2 2 J—k). As a result, for each value of 
k E {—2K, 0) U (0, 2K) there is only one negative mode. For these two bound states all the 6 field 
amplitudes in Eqs. dA. 141) — (IA. 191) are coupled, but in the limit q — > only 2 non-trivial amplitudes 
remain and the solutions reduce to that for the two decoupled equations (IA.22I) . 

Figj5] shows the dispersion relations uj^(k) and u%Jk) for several values of q. For q = they 
are given by the formula 

u 2 = (k ± Kf - K 2 (43) 

which will be explained shortly. For q > this formula is no longer exact, but it approximates 
reasonably well the values of uj 2 (k). (This formula was used to trace the q = 0.45 and q = 0.55 
curves in FiglU since the direct calculation of lo 2 (k), performed for q = 0.15, becomes rather 
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Figure 5: Left: dispersion relations lo 2 (k) for the two bound state solutions of Eqs. (IA.14l )-( |A.19b for the 
N = 1 twisted backgrounds with /? = 2. Right: The minimal value of ui 2 (k) corresponding to |k| = K. 
We see that the instability region bounded from below by the uj 2 (k) curve shrinks when q increases. 



time consuming for larger values of q.) The following features of this formula remain true also 
for q > 0. The negative mode exists only for k E (— 2K, 0) U (0, 2K). The minimal value of 
uj 2 is achieved for k = ±K and this value tends to zero as q — > 0, being bounded from below 
by — K 2 (see Fig©. We observe that the size of the instability region rapidly shrinks for q — > 1. 
The twisted vortices become therefore 'more stable' for large currents. This can be understood by 
noting that their energy decreases with growing X and for X — > oo approaches the lower bound 
27r|iV| Iil2ll and so that less room for instability is left for large X. 



Table I: The value of to 2 for k = K, N = 1. 



Q 


0.05 0.1 0.15 


0.2 


0.25 


0.3 


K 


0.3298 0.3294 0.3282 0.3259 


0.3225 


0.3181 0.3126 


u 2 


-0.108 -0.107 -0.105 -0.099 


-0.092 


-0.085 


0.076 


q 


0.35 0.4 0.45 0.5 


0.55 


0.6 


0.65 


K 


0.3059 0.2976 0.2876 0.2759 


0.2617 


0.255 


0.2259 


u 2 


-0.065 -0.055 -0.042 -0.031 


-0.0207 


-0.013 


-0.007 
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V. PHYSICAL MANIFESTATIONS OF INSTABILITY 



Summarizing the above analysis, despite their lower energy, the twisted vortices have essen- 
tially the same instabilities as the embedded ANO vortices. However, there is one important 
difference leading to important consequences: they do not have the uniform spreading instability 
mode independent of z analogues to the one obtained by setting k = in Eq.(fT6l). To under- 
stand how this comes about and why this is so important, let us explicitly reconstruct the fields 
corresponding to the unstable modes of the fundamental N = 1 background vortex. 



A. Reconstructing the perturbations 



Let tyr s \ with s = 1,2 denote the two linearly independent normalizedbound state solutions of 
the eigenvalue problem (|36l) for a given k and with m — 0, the corresponding eigenvalues being 
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. Each of them determines the 6 field amplitudes <pa \ Ua \ Q , xi m Eqs. dA. 141) — (1A. 191) and 



these, using (1341) . determine also ^\ Let us introduce vectors 



( d)?> K \ 



y 



V, 



a 



td,K I 
W ) 



( K' K \ 



(44) 



where \i = 3, 4 and where the components of the vectors Y^ )R and Z W)K are the m = amplitudes 
entering the mode decomposition (|27l) . If K describes a bound state solution of the perturbation 
equations, then it should be a linear combination of Xty since the latter can be regarded as the 
basis vectors in the space of bound state solutions. Similarly, if Z W:K describes a bound state 
solution of the second group of perturbation equations, then it should be a linear combination of 
the corresponding basis vectors. The latter, in view of the symmetry (|30l) . can be chosen to be 



X^ s y As a result, 



(45) 



where and bfy are independent real coefficients. Since the value of u is fixed for the bound 
states, the frequency assumes only two values, uj = ±VuA and one has similarly 



(»)' 



(46) 
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Inserting this to (1271) and applying (1241) gives the perturbations in the regular gauge, 

ft s 

S9 m = e iM^ iKz £ £ ^.hq^ + 



SA ? = EE + # } > A* = 3,4. (47) 

K S 

Here /ii s)± = i(0i s) ±^ s) ) and 

Qf s ) = ( A (s) exp{^ (s) t} + £f s) exp{-iu; (s) t}) exp{iKz} (48) 

with A* = al — %b K s and 5^ = c K s + id*. Regularity of the perturbations (1471) at p = is guaranteed 
by the boundary conditions (1381) . 



B. Instabilities of the embedded ANO vortex 

Let us consider again the ANO limit. We know that the bound states are then described by so- 
lutions of the two decoupled equations (|A.22I) . while Eqs. (IA.21l) do not have bound state solutions 
for N = 1. As a result, one has h^^ 1 = X2 = £m = an d> f° r s = 1, 

/4 1)+ = MP), h ( i ] - = 0, u 2 {1) = k 2 + - K\ (49) 

while for s = 2 

/4 2)+ = 0, 4 2) - = tt (p), c^fci-tf 2 , (50) 

where ^o(p) is me same as in Eq.(fT3l and k± = k ± K . This explains, in particular, the formula 
(1431 . The only non-zero perturbation in Eqs.(l47l) is then 

or explicitly 

5$ 2 = ^o(p) E ex P VV^*} + B W exp *} ex P{*M} 

K 

(Af 2) )* exp {-zy^t} + (5 ( K 2) )* exp {i^/o^t}] exp{-^}) . (52) 
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This expression contains exponentially growing in time terms if k± — K 2 = k(k ± 2K) < 0, 
which condition is satisfied for k E (-2K, 0) U (0, 2K). Since these terms are proportional to 
exp(±ik±z), the instability can be viewed as a superposition of standing waves of wavelength 
A = 2ir/k± whose amplitude grows in time. The minimal wavelength is A m i n = 2ix / K and this 
suggests that the instability can be removed by imposing the periodic boundary conditions along 
z-axis with the period L < A min . However, this will not remove the particular unstable modes with 
k± — independent on z since it can be considered as periodic with any period. Setting in (l52l) 
k + = or k_ = shows that these modes are proportional to each other and so in fact there is 
only one such homogeneous mode, 

5$ 2 = Ay (p)exp(Kt), (53) 

with A E C. This mode describes a uniform spreading of the vortex in the x, y plane. 

Terms with \k\ > 2K in (l52l) describe waves travelling towards positive and negative values of 
z with positive or negative frequency. These modes correspond to stationary deformations of the 
embedded ANO vortex by the twist. In particular, setting \k±\ — K and so k = 0, ±2K gives zero 
modes corresponding to static deformations of the vortex by the twist/current, 

<5$ 2 = tf (p)(4i e lKz + A 2 e- lKz ) J (54) 

two independent modes here corresponding to two possible directions of the current. 

It is worth noting that the elementary waves in (l52l) decouple one from another. In particular, 
denoting k = k + one can adjust the parameters A 1 ?, and Bfa such that the result will agree with 
Eq.©. 

C. Generic case 

Let us now consider perturbations of generic twisted vortices given by (|47l) . As in the ANO 
limit, there are two linearly independent bound state solutions of Eqs. dA. 141) — (1A. 191) with eigen- 
values ct>n\(«) and u;^ (k), where uj^Jk) < for —2K < k < and u;^ (ft) < if < k < 2K 
(see Figj5]). The perturbations are therefore given by Eqs.(l47l) where the sums contain growing in 
time terms for k E (-2K, 0) U (0, 2K). (When k — > 0, ±2K the solutions approach zero modes 
analogues to (l54l) and describing static on-shell deformations inside the family of twisted vortex 
solutions. In the generic case such zero modes contain a term linear in z, as can be seen by differ- 
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entiating the background fields © with respect to K, and so the solutions approach these modes 
only pointwise and non-uniformly in z.) 

The novel feature as compared to the ANO case is that now Eqs.(l47l) do not contain modes 
independent of z. This is related to the fact that Eqs. dA. 14h — (1A. 191) no longer decouple into inde- 

(V)4- (s) (s) 

pendent subsystems and all of the 6 radial field amplitudes h a , %g , £4 316 now non- vanishing. 



Terms with different 2-dependence in (1471) now turn out to be coupled to each other, which is the 
consequence of the fact that the background fields © depend explicitly on z. For example, it is 
no longer possible to have, as in Eq.(|49l), a solution with 7^ but with h 2 ~ = 0, and so the 
analog of the solution (1531) will read 

5$ 2 = A Vo(p)e Kt + B{1, t, p)e ±2lKz , 

where B(X,t,p) 7^ if X 7^ 0. Since all the unstable modes depend on z, the homogeneous 
negative mode is therefore absent for twisted strings 

Let us consider the effect of the perturbations on the gauge invariant vortex current. With Eq.© 
one obtains 

SJ Z = 9l(i6$* 2 d z $ 2 + i$* a d z 5§2 + 2A X $* 2 6$ 2 + f*6A x ) (55) 
and using Eqs.(|6l).(l47T) gives the growing part of this expression in the form 

2K 

SJ z = f 2 (p) Yl ^(A(p)cosM + S K (p)sinM), (56) 

k=-2K 



where uj = \ for k < and uj = y \ for k > 0. A similar expression can be obtained for 

5Jq. Every unstable mode with k E (-2K, 2K) therefore produces ripples with the wavelength 
A = 27r/|«| on the homogeneous distribution of the background current density. As the amplitude 
of these ripples grows in time, the current deforms more and more thus tending to evolve into a 
non-uniform 'sausage like' structure characterized by zones of charge accumulation. 

Since the instability exists only for a finite range of k, there is a minimal instability wavelength 
corresponding to the maximal value of k, 

IX 

Amin T r • (57) 
A 

As a result, finite vortex pieces obtained by imposing periodic boundary conditions on the infinite 
vortex will have no room for inhomogeneous instabilities if their length L is less than A max . Next, 
as we already know, twisted vortices do not have homogeneous negative modes, which means 
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Figure 6: Schematic view of the perturbed vortex current. 



modes periodic with any period (for example, the k = mode in (|56l) is homogeneous but not 
negative). It follows then that short vortex pieces do not have negative modes at all. They are 
therefore stable. 

At first view the condition L < ir/K is impossible, since the background vortex configuration 
© contains exp(iKz) and so L should be an integer multiple of 2n/K. However, if the semilocal 
model is viewed as a subsector of the local SU(2) x U(l) theory in the limit where the SU(2) 
gauge field decouples, it is possible to use the local gauge transformations (I2TT) not merely as an 
intermediate technical tool as was done above, but in order to pass to and stay in a physically 
admissible gauge where the vortex fields do not depend on z. For example, applying to © the 
gauge transformation (|2TI) generated by 



This gauge is globally regular and contains no explicit z dependence. One can therefore always 
work in this gauge and identify z with any period. 

To recapitulate, short pieces of twisted N — 1 vortices are stable, but they become unstable 
with respect to fragmentation if their length exceeds a certain critical value. It is worth mention- 
ing an interesting analogy with the well known Plateau-Rayleigh instability of a fluid cylinder in 
hydrodynamics (see e.g. J3]). The cylinder is stable if only it is short enough, but when it gets 
longer the surface tension tends to break it up into disjoint pieces. This can easily be seen when 
water comes out of a kitchen tap: while close to the tap the water jet is homogeneous, far enough 




gives (setting for simplicity = M = 0) 
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from it the surface ripples appear. This hydrodynamical analogy seems to be not accidental, since 
the twisted superconducting vortex can effectively be described as a superposition of two charged 
fluids flowing in the opposite directions J12Q. Another interesting analogy that could be mentioned 
comes from a completely different domain: the Gregory-Laflamme instability of black strings in 
the gravity theory in higher dimensions (2(X\ . Black strings become unstable with respect to inho- 
mogeneous perturbations if their length exceeds a certain critical value, which can be explained 
by the tendency of the event horizon to minimize its area thus reducing the black string entropy. 

It is worth mentioning that in both of these two examples, and also in our case, the unstable 
modes exist only in the axially symmetric m = sector. 



VI. CONCLUDING REMARKS 



Summarizing what has been said above, we considered generic field fluctuations around the 
twisted current carrying vortices in the SU(2) g i ob ai x U(l)i oca i semilocal field model, or equally in 
the electroweak theory in the limit where the weak mixing angle is n/2. By studying the negative 
modes we have concluded that the twisted vortices exhibit essentially the same instabilities as the 
embedded ANO vortices, but with one important exception: they do not have the uniform spread- 
ing instability. All the negative modes of the fundamental N — 1 twisted vortex are non-uniform, 
with the wavelength bounded from below by ir/\K\, where K is the vortex twist parameter. Since 
the uniform instability is absent, it follows that short enough vortex pieces should be stable, since 
they have no room to accommodate the non-uniform longwave instabilities. Longer vortex pieces 
will be unstable and the instability will lead to their fragmentation into non-uniform objects char- 
acterized by zones of charge accumulation and by a non-uniform current density. However, to 
fully describe this process requires going beyond the linearized approximation. 

It is interesting to note that increasing the vortex current X decreases the twist K (see FigH]) 
and so that the minimal instability wavelength k/\K\ increases, which makes stable longer and 
longer vortex pieces. In addition, the characteristic time of the instability growth, 2tt/\uj\, also 
increases, since < \K\ (see Figj5]). As a result, increasing the current removes more and more 
of negative modes thus making the vortex 'more and more stable' . 

The non-uniform character of perturbations suggests that if there is a stable configuration to 
which the perturbed vortex could finally relax, then this configuration should be inhomogeneous. 
The current should be homogeneous in such a state, though. Specifically, integrating the current 
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conservation condition = over a 3-volume sandwiched between z = Z\ and z = z 2 planes 
gives 

d ^=X{z l )-X{z 2 ), 

where Q is the total charge in the volume and T(z\) is the total current through the z — Z\ 
plane, similarly for T{z,2). The condition T{z\) ^ 1{z%) therefore requires that ^ 7^ 0, and 
so inhomogeneities of the current can only appear during the dynamical phase when the time 
derivatives are non-zero. However, if the perturbed vortex finally ends up in an equilibrium state 
which is stationary, then one will have T{z\) = X(z 2 ) and so the current will be homogeneous. At 
the same time, as we already know, the vortex itself cannot be homogeneous in a stable state (since 
the homogeneous vortex is unstable). We therefore conclude that if a stable stationary state of a 
current carrying vortex exists, then in this state the current will be constant along the vortex, but 
the vortex configuration itself will depend non-trivially on z and perhaps also on (p. An interesting 
problem would be to look for such a stable state via minimizing the energy by keeping the current 
and also total charge fixed. 

Another possibility to have stable objects would be to take short vortex pieces and to close them 
to make small loops. If exist as stationary field theory solutions, they might be stable. The issue of 
constructing such solutions is actually quite important, since stationary vortex loops stabilized by 
centrifugal force (they are called vortons |2l|]) have been extensively discussed for about last 20 
years (see [2J] for a review). However, to the best of our knowledge, such objects have never been 
constructed by resolving the field equations of the underlying field theory, excepting two examples 
obtained in a rigid and not local gauge field theory 12211 . 

It is finally worth noting that the above results are likely to be generalizable in the context of 
the full electroweak theory for arbitrary values of the weak mixing angle, since the analogs of the 
twisted vortices are known to exist in this case \ 2^\ . It seems plausible that short pieces of these 
superconducting electroweak vortices could also be stable, and also perhaps small vortex loops. 
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Appendix 



We derive here the full system of the perturbation equations. The starting point is the mode 
decomposition (|27T) . inserting which to the fluctuation equations (|25l) and separating the t, z, ip 
variables gives, for each set of values of cu, k, m, a system of 16 ODE's for the 16 radial field 
amplitudes in (1271) . These equations split into two independent subsystems of 8 equations each 
which have exactly the same structure, upon the identification (l30l) . Equations of the first group 
involve the amplitudes <f> a , u a , xa, Xs> Xa- 
N 2 ( ai - lf + m 2 



E 1 



-V, 



P 



+ P(3ft + K - 1) + (K 2 - n 2 )a 2 2 + k 2 -lu 



: 2 if L j 2 o,+ 1 ( ( ihs Kh)u, + Nm{l 2 ai) ) //, - 2<>r, 2 /jsJ + 2Ku,h 



P 



2Nf 1 (a 1 -l) 
P 2 



U = o, 



(A.l) 



-V + + N2(1 a f + m \ PUl + fl ~ 1) + (K 2 - n 2 )aj + n 2 -u 2 



P 



+2 ({Sk* - Kn)a 2 + Nm{X Ql) ) fa + w/i 6 + /1X2 + ^ + 



X'2 



6 r S4 = 0, 

p 2 



(A.2) 



-V + + (iVai A f + m2 + (3{fl + 3/| - 1) + (7T 2 - n 2 )(a 2 - I) 2 + k 2 cu 2 
P 

+2f3hf 2 fa + 2 f (Qu; - K«)(a 2 - 1) + m ( M ~ Na ^ v% + 2 fi/ 2 (i - a 2 ) & 



+ 2^ 2 (a 2 -l)e3 + 2/2(iVa r M) ^0, 



(A.3) 



4 •- 



A + (iVfll M / +m2 + - fi 2 )(a 2 - I) 2 + + f 2 -l) + n 2 -u 2 



P 



V-2 



+2 (Qui - KK)(a 2 - 1) + 



m(M - Nai 
? 



' 2 + a;Mi + /2X2+(2/ 2 + y jxs 



-«/ 2 ^ 3 



m/ 2 



6 = 0, 



(A.4) 



m 



+ — + 2(/ 2 + / 2 )-^ 
P 



e 3 + 4K/ 1 a 2 1 + 4ftr/ 2 (a 2 -l)</. 2 



-2k/i Pi - 2k f 2 i/ 2 + wfi:(i + dx2 + - - 6 = 0, 

p J p 2 



(A.5) 
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(A.6) 



E 7 := -V + + — + 2{fl + f 2 ) + K 2 ^ 1 +AQf 1 a 2 <p 1 + AQf 2 (a 2 -l)<p 2 



-2w/i Vi - 2ujf 2 u 2 + u [x' 2 H 6 5- ^4 = 0, 

VP/ P 




(A.7) 



E 8 




- + 2(/ 2 + / 2 2 ) + - cu 2 )x2 + 2(/>i - M) + 2(/^2 - M) + < 



(A.8) 



Here "D-t := ^ ± the functions /i,/ 2 ,ai,a 2 and the constants Q,K,N,M refer to the 
background twisted vortex solution. 

Not all of these equations are independent. Using the background field equations © one can 
check that the following identity holds, 



and so one of the equations is redundant. 

A good consistency check for the equations is provided by the translational modes. If 
A^x 01 ), <fr(x a ) is the background vortex solution written in the gauge ©, then A^(x a + 
X a ),$>(x a + X a ), where X a is a constant vector, will also be a solution. This implies that 
the Lie derivatives of A^(x a ), &(x a ) along X a , 



should fulfill the linearized field equations. If X — the vector generating displacements of 
the vortex in the x direction, then calculating the Lie derivatives (|A.10I) and transforming them 
according to ((24)) gives the perturbations in the form (HT)) with m = l,u = k = and 



E' s H E e — kE 5 + uE 7 - 2f x E 2 - 2f 2 E A = 0, 



P P 



(A.9) 



SAp = X a O a A^ + A^X", 5$ = X a d a <$> 



(A.10) 




-a 2 ). 
P 



(A.ll) 
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Using the background field equations © one can check that these fulfill Eqs. dA. 1 1 )— (1A.8I) . Choos- 



ing similarly X = gives uj = n = m = and 



b a = 0, v x = 0, v 2 = / 2 , €i = X2 = €3 = U = 0, 



(A.12) 



which also solves the equations. 

Using again Eqs.© one can check that Eqs. dA. 1 1 ) — (IA.8I) are invariant under the gauge transfor- 
mations (T29l) . 



*>a, "a -> Va + faJ, 6 ~> £l+^7> X2 ~» X2+7', 6 ~> &+«7, & ~> ^4 + ^7, (A.13) 



where 7 = 7(p). In order to fix this gauge symmetry we impose the temporal gauge condition 
5A = by setting £1 = 0. We also specialize to the case of purely magnetic backgrounds, 
fi = 0. Eq. (IA.7l) can then be resolved with respect to £ 3 (see Eq.(l34l). At the same time one 
can exclude from consideration Eq. (IA.5l) containing derivatives of £ 3 , since we know that one of 
the equations in the system is redundant. Inserting then (1341) into the remaining equations (IA.1I) - 
(IA.4I) . (|A.6I) . (IA.8I) gives the 6 independent second order equations for 6 field amplitudes a , v a , X2, 
£4, 

N 2 ( ai - lf + m 2 



-V + + 



+ /5(3A i + / 2 2 -l) + K 2 a^ + ^-c 



, ( Nm(l - 01) u 2Ka 2 ff \ 4Ka 2 f\f 2 1Kj x a 2 ( X i 

+ 2 KKa 2 vx u 2 H [ x 2 H 

V p K J K 



K 



+ ^ - i) 



Kma-: 



£4 = 



(A. 14) 



-p+ + ^ a . l)2 + m2 + KA 2 + / 2 2 - 1) + ^ + VI + k 2 -^ 2 



| / iVm(l- ai ) 
V P 2 



#«a 2 J 0i + 2/ 1 / 2 z/ 2 + 2/( X 2 = O, 



(A.15) 



- (A^ai-M) 2 + m 2 2 2 2 2 2 2 
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+ 2/3/i/2. 



4if(l - agjA^ n f m(M-N ai ) T ^ 2K(l-a 2 )f 2 2 \ 
H ^1 + 2 ( 5 h ak(1 - a 2 ) H ] z/ 2 



+ 



2K/ 2 (a 2 -r 



X 2 



P z 

*A , 2/ 2 



iVai - M + 



Km(l — a 2 ) 



£4 = 



(A.16) 
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-£>+ + ( (M Na f + m2 + /3(/ 2 + /f - 1) + K\l - a 2 ) 2 + 2ft + k 2 -uj 2 

p ~ . 

-I m( < M - Na ^ + - a 2 )) 2 + v x + 2f> X 2 = 0, 



v 2 

(A.17) 
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D_ + — + 2/ 2 + 2/ 2 + - 



2 777, 

£ 4 +4iV/i( ai - l)0i + 4/ 2 (iVai - M)0 2 X 2 = 
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(A. 18) 



m 2 + 1 
7^ 



+ 2/ 2 + 2/| + k 2 - cu 2 



X2 + 4/>! + 4/> 2 



2m 

7 



o. 



(A.19) 



As explained in the main text, there is no gauge freedom left in this system. These equations are 
invariant under 

m -> -m, k -> ^ 2 -> u; 2 , a -> a , v a -> -i/ a , £ 4 -> £ 4 , X2 -> -X2 • (A.20) 

Taking suitable linear combinations of the 6 field amplitudes a , z/ a , % 2 , £ 4 and their first derivatives 
one can build a 6-component vector $ using which the 6 equations (IA. 141) — (1A. 1 91) can be rewritten 
in the form of a Schrodinger-type eigenvalue problem d36l) . 

In the ANO limit Eqs. dA.14b - dA.19b split into two independent subsystems. The first subsys- 
tem is obtained by setting in Eqs. dA. 14I) . (|A. 1 5I) . (|A. 1 8I) , (|A. 191) f 2 = a 2 = 0: 

2N(1 - ai) 



- A^ 2 («i - l) 2 + m 2 „, - . , 9 
-X>+ + 1 ,/ + /?(3/i - 1) + k 2 - cu 2 
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iV 2 (l-ai) 2 + m 2 



>i + 



(mz/i - /i£ 4 ) = 



+ /3(/ 2 -1) + 2/ 2 + k 2 -^ 



2iVm(l - 01) , 
"i+ — -4>i + 2f[ X 2 = 



m 



-V_ + ^-+ 2 ft + K 2 -LU 2 

p 

m 2 + 1 



9777 

e 4 + 4iV/ 1 (a 1 -l)0 1 X 2 = 0. 

P 



-V + + 



P' 



+ 2ft + k z -uj' 



. 2m 

X2 + v x =- £ 

p 2 



0. 



(A.21) 



Since these equations do not contain perturbations of the second component of the Higgs field, 
they actually describe the dynamics of the perturbed ANO vortex within the original ANO model. 
Eqs. dA.16b . dA.17b in the ANO limit reduce to two decoupled equations for the perturbations of 



the second component of the Higgs field h 2 



± v 2 



V + + ^ 2 + } + f3(f 2 -l) + (K± k) 2 

p 2 



to 



(A.22) 
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